Environmental  Laboratory  erdc/el  tr-05-5 


US  Army  Corps 
of  Engineers  ® 

Engineer  Research  and 
Development  Center 


Electromagnetic  Power  Attenuation  in  Soils 

John  O.  Curtis  August  2005 


Approved  for  public  release;  distribution  is  unlimited. 


ERDC/EL  TR-05-5 
August  2005 


Electromagnetic  Power  Attenuation  in  Soils 

John  O.  Curtis 
Environmental  Laboratory 

U.S.  Army  Engineer  Research  and  Development  Center 
3909  Halls  Ferry  Road 
Vicksburg,  MS  39180-6199 


Final  report 

Approved  for  public  release;  distribution  is  unlimited 


Prepared  for  U.S.  Army  Corps  of  Engineers 
Washington,  DC  20314-1000 


ABSTRACT:  This  study  explores  the  possibility  of  predicting  the  attenuation  of  electromagnetic  energy  in  soils 
as  a  function  of  soil  type,  soil  moisture  content,  and  the  frequency  of  the  sensor  system  being  used  to  probe  the  soils. 
Of  primary  concern  is  the  issue  of  safety;  one  does  not  want  to  risk  an  explosion  when  using  high-power  electro¬ 
magnetic  sensors  to  search  for  buried  unexploded  ordnance.  While  explicit  predictive  models  could  not  be  found 
because  of  the  unpredictable  effects  of  soil  chemistry,  simple  statistical  analyses  did  provide  conservative  guidance 
at  fixed  frequencies  and  selected  soil  moisture  contents. 
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Introduction 


Problem 

This  study  investigates  the  suitability  of  simple  mathematieal  models  for 
predieting  the  attenuation  of  eleetromagnetie  (EM)  energy  in  soils.  The 
motivation  for  this  effort  is  a  safety  issue.  Explosives  and  ordnanee  (EOD) 
eleanup  erews  may  need  to  utilize  eleetromagnetie  sensors  to  deteet  and  identify 
unexploded  ordnanee  (UXO)  buried  in  various  types  of  soils.  They  must  be 
reasonably  eertain  that  their  instruments  will  not  trigger  a  premature  detonation 
of  the  UXO. 

Experienee  has  shown  (see  diseussion  in  the  next  seetion)  that  EM  power 
attenuation  in  soils  is  a  strong  funetion  of  the  ffequeney  of  the  waves,  the  amount 
of  moisture  in  the  soil,  and  the  type  of  soil  (soil  ehemistry).  Soil  density  and  soil 
temperature  are  less  important  faetors.  Therefore,  the  goal  of  this  effort  was  to 
generate  mathematieal  models  to  prediet  EM  attenuation  as  a  funetion  of  soil 
type,  soil  moisture,  and  frequeney  of  the  sensor  system.  If  EOD  erews  (or 
seientists  and  engineers  who  are  testing  new  equipment)  are  able  to  elassify  field 
site  soils  as  well  as  estimate  the  volumetrie  moisture  eontent  of  the  soils,  then  the 
models  pursued  in  this  study  ean  be  utilized  to  estimate  power  attenuation  in 
those  soils  and  identify  eonditions  under  whieh  the  erews  might  faee  risks  in 
working  with  the  equipment. 


Theory 

A  plane  eleetromagnetie  wave  traveling  through  soil  whose  eleetrieal 
behavior  ean  be  eharaeterized  by  a  frequeney-dependent  complex  dielectric 
permittivity  will  experience  a  loss  of  power  that  can  be  calculated  by  the 
following  formula  (see  Appendix  A): 


power  loss  ( 
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where 


s' ,  e"  =  real  and  imaginary  parts  of  the  relative  eomplex  dieleetrie 
permittivity,  respeetively 

m  =  angular  ffequeney  of  the  wave  ( =2itf  ,f\%  the  ffequeney  in  eyeles 
per  seeond,  or  Hz  ) 

c  =  speed  of  light  (3x10*  m/s). 

The  ratio  of  the  imaginary  part  of  the  dieleetrie  permittivity  to  the  real  part  is 
often  referred  to  as  the  loss  tangent. 

If,  instead  of  eharaeterizing  a  soil  by  its  eomplex  dieleetrie  permittivity,  it  is 
deseribed  by  a  real  dieleetrie  permittivity  and  an  eleetrieal  eonduetivity,  then 
Equation  1  eould  be  written  as: 


power  loss 
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where 


o  =  eleetrieal  eonduetivity  (mho/m,  or  seimens/m) 

E()  =  eleetrieal  permittivity  of  free  spaee  (8.85x10’'^  farads/m) 

Equation  2  ean  be  used  to  generate  power  attenuation  predietions  for  given 
values  of  the  dieleetrie  eonstant  and  eleetrieal  eonduetivity.  The  results  are 
shown  in  Figures  1  and  2  on  the  following  pages.  In  Figure  1,  power  attenuation 
has  been  plotted  against  operating  ffequeney  for  a  eonduetivity  value  of 
1 0  milliSiemens  per  meter  (mS/m),  whieh  is  typieal  of  low  eonduetivity  soils, 
and  various  values  of  the  dieleetrie  eonstant.  Figure  2  is  the  same  ehart  for  a 
eonduetivity  value  of  50  mS/m,  whieh  eould  be  representative  of  high 
eonduetivity  soils.  These  eharts  ean  only  be  used  if  one  knows  the  values  of 
dieleetrie  eonstant  and  eleetrieal  eonduetivity.  The  problem  is  that  the  eleetrieal 
properties  of  soils  are  a  funetion  of  ffequeney,  moisture  eontent  of  the  soil,  the 
ehemistry  of  the  soil  that  determines  how  many  ions  ean  be  produeed  when 
moisture  levels  ehange,  the  temperature  of  the  soil/water/air  matrix,  and,  in  some 
very  eomplex  manner,  the  strueture  of  the  matrix,  whieh  is  related  to  its  dry 
density,  its  porosity,  and  the  type  of  minerals  it  eontains. 

Another  interesting  point  to  make  from  the  theoretieal  predietions  shown  in 
Figures  1  and  2  is  that,  eontrary  to  intuitive  understanding,  eleetrieal  eonduetivity 
alone  does  not  eontrol  EM  power  loss  in  soils.  Clearly,  all  other  things  being 
equal,  power  attenuation  inereases  with  inereasing  eonduetivity.  However,  as 
these  figures  show,  at  some  frequeneies,  a  low  dieleetrie  permittivity,  low 
eonduetivity  soil  ean  experienee  as  mueh  power  attenuation  as  a  high 
permittivity,  high  eonduetivity  soil.  While  it  may  seem  far-fetehed,  a  possible 


2 


Chapter  1  Introduction 


Conductivity  =  10  mS/m 


Frequency  -  MHz 


Figure  1. 
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explanation  is  that  the  low  permittivity  soil  allows  the  eharged  partieles  more 
opportunity  to  be  exposed  to  eleetrieal  loss  meehanisms  (by  not  allowing  them  to 
be  trapped  as  stored  energy),  thus  effeetively  produeing  more  eleetrieal  loss. 
What  this  means  is  that  one  must  be  very  eautious  in  utilizing  “rules-of-thumb” 
for  power  attenuation  in  soils  based  on  field  measurements  of  effeetive 
eonduetivity. 


Previous  Soil  Property  Models 

Clearly,  the  problem  of  predieting  EM  attenuation  in  soils  reduees  to 
predieting  the  eleetrieal  properties  of  the  soils,  whether  they  be  expressed  as  a 
eomplex  dieleetrie  permittivity  or  a  real,  relative  dieleetrie  permittivity  and  an 
eleetrieal  eonduetivity.  The  literature  is  replete  with  attempts  to  model  soil 
eleetrieal  properties  as  funetions  of  other  soil  parameters  sueh  as  moisture 
eontent,  soil  texture,  and  resistivity  of  pore  water,  ete.  Unfortunately,  at  best, 
these  models  are  site  speeifie,  and  one  eannot  avoid  scooping  up  samples  in  the 
field  that  must  be  analyzed  in  a  laboratory. 


Archie’s  Law 

As  an  example,  consider  Archie’s  Law  (Archie  1942)  for  calculating  the 
effective  electrical  resistivity  of  a  soil: 

=  acp-s-'^p^  (3) 

where 

(p  =  the  fractional  pore  volume  (porosity) 
s  =  the  fraction  of  pores  containing  water 
=  the  resistivity  of  the  pore  water 
a,  n,  and  m  =  constants 

This  formula  gives  the  soil  conductivity  (the  inverse  of  the  resistivity),  but  only 
after  conducting  laboratory  measurements  of  soil  porosity  and  pore  water 
resistivity  and  doing  whatever  is  necessary  to  estimate  the  other  parameters. 
Clearly,  direct  measurements  of  soil  attenuation  in  a  laboratory  are  simpler. 


Army  study 

Another  early  model  of  soil  conductivity  came  from  an  extensive  field 
measurement  program  conducted  by  the  U.S.  Army  (Hulse  et  al.  1972).  An 
equation  for  conductivity,  a,  was  developed  for  regional  soils: 

cr  =  1  +  8.4(w  ;  -  0.  Ic  (4) 
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where 


c  =  the  elastie  fraetion  by  pereent  weight  (assumed  to  represent  the 
fraetion  of  soil  eontaining  mineral  and  roek  fragments) 

m:k  =  weight  ratio  of  montmorillonite  to  kaolinite  elays  in  the  soil. 

Using  this  model  requires  extensive  soil  ehemistry  experiments  to  determine 
elay  mineral  ratios,  as  well  as  whatever  teehniques  are  needed  to  estimate  the 
elastie  fraetion. 


University  of  Kansas  studies 

Possibly  reeognizing  that  sueh  models  are  too  esoterie  for  many  field 
engineers  and  seientists,  a  group  at  the  University  of  Kansas,  and  later  at  the 
University  of  Miehigan,  attempted  to  model  soil  eleetrieal  properties  as  a 
funetion  of  soil  physieal  properties  (Hallikainen  et  al.  1985).  The  following 
formulae  were  proposed  for  values  of  the  eomplex  dieleetrie  eonstant: 

s'  =  {g'q  +  a'^S  +  a2C)+  {b'^  +  b'^S  +  b'2C)m^  +  (c^  +  CjiS  +  (5) 

s"  =  [a'l  +  a'IS  +  a"c)  +  [b'l  +  b'^S  +  +  (c"  +  c'^S  +  c" c)w'  (6) 

where 

S  =  the  weight  fraetion  of  sand  partieles  in  the  soil  (diameter  between 
0.05  mm  and  2.0  mm) 

C  =  the  weight  fraetion  of  elay  partieles  in  the  soil  (diameter  less  than 
0.002  mm) 

niy  =  the  volumetrie  moisture  eontent  of  the  soil 

and  the  eoeffieients  are  determined  by  fitting  laboratory  data  eolleeted  over  the 
range  of  frequeneies  being  studied.  Eaeh  frequeney  has  assoeiated  with  it  a 
different  set  of  eoeffieients. 

These  equations  are  exaetly  what  is  desired  for  this  study:  seleet  an 
operating  frequeney;  make  a  site  measurement  of  moisture  eontent;  perform  a 
simple  onsite  soil  elassifieation  to  determine  sand  and  elay  fraetions;  and 
ealeulate  the  eomplex  dieleetrie  eonstant.  From  Equation  1 ,  one  ean  then  prediet 
EM  attenuation.  Unfortunately,  the  eoeffieients  are  soil-dependent.  In  faet,  it  is 
highly  probable  that  two  soils  with  similar  texture  and  moisture  values  will 
produee  vastly  different  dieleetrie  eonstant  values  beeause  of  differing 
ehemistries. 


Electrical  network  modeling 

Another  approaeh  to  modeling  soil  eleetrieal  properties  that  aeeounts  for 
their  dispersive  (frequeney-dependent)  nature  is  that  of  using  eleetrieal  networks 
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of  resistors  and  capacitors  whose  parameters  are  chosen  to  fit  specific  laboratory 
data  (Smith  and  Arulanandan  1981).  For  example,  the  apparent  conductivity  of 
the  soil  can  be  modeled: 


d{\-d)S 


KK  k  k  2  2 

\-d  d 


\-d 


J 


+  bk^  +  ck. 


(7) 


in  which 


S  = 


- i - 1 - L- 

\-d  d 


2  2 

+  G)  s  i 


f 

S, 

^  _| _ L 

yl-t/  d  J 


a,  b,  c,  and  d  are  geometrical  weighting  parameters  associated  with  electrical  path 
lengths  through  a  soil  sample  and  relative  amounts  of  soil  solids  and  soil 
solutions,  while  Sr  and  kr  are  the  relative  dielectric  constant  and  electrical 
conductivity  of  the  soil  fraction,  and  Ss  and  ks  are  the  relative  dielectric  constant 
and  electrical  conductivity  of  the  soil  solution. 


Once  again,  a  model  such  as  Equation  7  is  not  an  improvement  over  simply 
making  permittivity  and  conductivity  measurements  directly.  The  weighting 
parameters  are  found  from  real  conductivity  data,  anyway.  Furthermore,  every 
different  soil  has  a  different  set  of  weighting  parameters. 


Other  studies 

There  are  many  other  references  in  the  literature  to  soil  electrical  property 
measurements  and  modeling  (Flipp  1974;  Floekstra  and  Delaney  1974;  Campbell 
1990;  Marinelli  et  al.  1999;  Miller  et  al.  2002;  Cooke  and  Gladwin  2004).  As 
with  those  above,  they  all  fall  short  of  the  requirements  for  this  study  in  that  none 
of  them  provide  a  universal  model  for  all  soils,  or  even  broad  classes  of  soils. 

The  earlier  studies  always  dealt  with,  at  most,  a  few  different  soil  samples. 


Scope  of  Investigation 

The  hope  for  this  investigation  was  that  a  statistical  analysis  of  a  large 
amount  of  frequency-dependent  laboratory  data  on  soils  of  varying  texture, 
moisture  content,  bulk  density,  and  possibly  temperature  might  result  in  an 
acceptable  model  for  predicting  EM  attenuation  from  simple  measurements  that 
can  be  conducted  easily  in  the  field.  The  chapters  that  follow  document  this 
effort. 

Chapter  2  begins  with  a  brief  description  of  how  electrical  property  data  have 
been  collected  at  ERDC  and  goes  on  to  display  those  data  at  a  single  frequency 
and  a  single  temperature.  The  two  soil  classification  schemes  used  in  this  effort 
are  also  briefly  discussed.  The  chapter  concludes  with  a  visual  representation  of 
both  soil  classifications. 
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Chapter  3  steps  through  the  analysis  of  the  attenuation  data  at  eaeh  of  five 
seleeted  ffequeneies.  Arguments  are  made  for  whether  or  not  eaeh  of  the  soil 
physieal  properties  play  a  signifieant  role  in  predieting  EM  attenuation. 

Conelusions  are  summarized  in  Chapter  4,  and  reeommendations  are  made  as 
to  how  to  provide  a  faetor  of  safety  to  equipment  users  in  the  field. 

An  appendix  is  provided  that  ineludes  the  development  of  the  theoretieal 
relationships  among  the  eleetrieal  properties  and  various  physieal  properties  of 
the  soils  that  were  tested  at  ERDC. 
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ERDC  Soil  Data 


Since  about  1990,  the  author  has  measured  the  complex  dielectric  constant  of 
hundreds  of  soil  samples.  Many  of  these  samples  were  taken  from  military 
installations  in  the  United  States  as  well  as  from  locations  in  Europe  and  the 
Middle  East.  A  number  of  samples  of  “prepared”  soils  from  the  U.S.  Army 
Corps  of  Engineers  Geotechnical  Eaboratory  in  Vicksburg,  MS,  as  well  as 
several  benchmark  soils  from  the  National  Soils  Survey  Center  in  Eincoln,  NE, 
were  also  studied. 

Fine-grained  materials  were  packed  into  coaxial  sample  holders  whose 
temperature  could  be  controlled.  Normally,  the  sample  holders  were  packed  with 
the  soil  as  it  was  drawn  from  the  container  in  which  it  was  shipped,  and  a 
measurement  of  the  complex  dielectric  properties  of  that  sample  were  made  over 
a  selected  band  of  frequencies.  If  the  soil  was  noticeably  damp,  it  might  be  air- 
dried  while  still  in  the  container,  after  which  its  properties  were  measured  again. 
Finally,  distilled,  deionized  water  might  be  added,  and  the  sample  allowed  to 
reach  an  equilibrium  condition  before  another  measurement  was  made.  In 
summary,  each  physical  sample  would  be  interrogated  at  at  least  three  different 
moisture  contents  (Curtis  2001b). 

While  each  measurement  set  covered  different  specific  frequency  limits, 
most  of  the  data  cover  overlapping  frequency  bands.  For  this  reason,  and  given 
the  measurement  technique  described  above,  the  author  believes  these 
measurements  represent  a  statistically  significant  data  set  whose  independent 
variables  include  soil  type,  frequency,  volumetric  moisture,  soil  bulk  density, 
and,  to  a  lesser  extent,  soil  sample  temperature.  In  this  study  it  was  hoped  that 
such  a  data  set  would  be  a  viable  source  of  EM  attenuation  prediction  models. 


Single-Frequency  Data 

As  an  example  of  data  collected  since  1 990,  the  following  charts  contain  soil 
electrical  property  data  plotted  against  volumetric  moisture  content  at  one 
frequency,  100  MHz,  and  at  one  temperature,  20  deg  C.  There  are  approximately 
1,080  data  points  on  each  chart,  representing  between  300  and  400  distinct  soil 
samples.  It  is  also  important  to  note  that  these  scatter  plots  distinguish  neither 
soil  classifications  nor  soil  bulk  density.  Therefore,  the  only  variable  on  this  first 
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set  of  charts  is  volumetric  moisture  content.  The  next  chapter  will  contain  a 
study  of  attenuation  data  as  a  function  of  the  other  test  and  soil  parameters. 

Laboratory  data  are  first  processed  to  produce  frequency-dependent  values 
of  the  complex  dielectric  constant  using  the  algorithm  described  in  Appendix  A. 
Other  electrical  parameters  are  then  calculated  from  the  complex  dielectric 
constant  values  using  equations  that  are  also  found  in  Appendix  A. 

Figures  3  and  4  contain  the  real  and  imaginary  parts  of  the  complex  dielectric 
constant,  respectively.  Clearly,  there  is  a  very  strong  relationship  between  the 
real  part  of  the  complex  dielectric  constant  (which  is  also  referred  to  as  the 
dielectric  permittivity)  and  volumetric  moisture  content.  Other  than  a  few 
“outlier”  data  points,  one  could  readily  accept  a  model  fit  to  these  data  by 
regression  that  would  pass  through  the  center  of  the  narrow  band  of  data.  In  fact, 
this  is  exactly  what  is  done  with  some  field  equipment  used  to  measure  soil 
moisture.  If  one  is  willing  to  accept  the  spread  of  data  in  Figure  3  as  the  error 
bar,  then  a  model  can  be  easily  generated  that  predicts  soil  moisture  from  a 
measurement  of  the  dielectric  permittivity  (Topp  et  al.  1980;  Curtis  2001a).  A 
value  of  the  permittivity  can  be  obtained  from  a  field  instrument  that  treats  the 
soil  at  the  end  of  the  open-ended  coaxial  probe  as  a  shunt  capacitance. 

The  five  outlier  data  points  on  Figure  3  that  fall  above  the  band  of 
concentrated  data  all  come  from  measurements  on  a  manufactured  soil  that  was  a 
mix  of  clean  sand  and  magnetite.  It  is  entirely  possible  that  the  high  values  of 
dielectric  permittivity  are  a  product  of  the  assumption  used  in  data  processing 
that  the  soil  is  nonmagnetic.  One  of  the  four  outlier  data  points  that  fall  below 
the  concentration  band  is  from  a  soil  that  is  a  natural  magnetic  mineral  bearing 
soil.  The  chemistry  of  the  other  three  points  is  unknown  other  than  two  of  them 
visually  classified  (Unified  Soil  Classification  System)  as  an  organic  silty  sand 
from  Alaska,  and  the  other  is  a  soil  sample  taken  from  a  2-meter  depth  at  the 
Aberdeen  Proving  Ground,  MD.  In  fact,  these  low-value  outliers  may  represent 
nothing  more  than  a  computational/experimental  error. 

While  a  strong  relationship  exists  between  dielectric  permittivity  and 
volumetric  moisture  content,  virtually  no  relationship  exists  between  the 
parameters  that  reflect  energy  losses  in  the  soil  and  moisture  content.  For 
example,  the  imaginary  component  of  the  complex  dielectric  constant  (Figure  4) 
plotted  against  volumetric  moisture  produces  a  true  “scatter”  plot.  Both  the  loss 
tangent  and  the  apparent  conductivity  (Figures  5  and  6,  respectively)  produce 
similar  results. 

Although  power  attenuation  and  phase  velocity  (Equations  A16  and  A 17, 
respectively)  are  functions  of  both  the  real  and  imaginary  parts  of  the  complex 
dielectric  constant.  Figure  7  shows  that  power  attenuation  cannot  be  predicted  by 
moisture  content  alone.  On  the  other  hand.  Figure  8  reveals  that  the  phase 
velocity  of  EM  waves  traveling  through  all  of  the  soils  tested  at  ERDC  is  highly 
correlated  with  volumetric  moisture  content.  This  latter  observation  is  supported, 
in  part,  by  noting  that  for  the  majority  of  loss  tangent  values  measured  at  ERDC, 
the  phase  velocity  will  be  inversely  proportional  to  the  square  root  of  the 


Chapter  2  ERDC  Soil  Data 


9 


Figure  3. 


Real  (dielectric 
constant)  vs.  volumetric 
moisture  content, 

100  MHz,  20  deg  C 


All  ERDC  Samples 
100  MHz  Data 


volumetric  moisture  -  % 


All  ERDC  Samples 
100  MHz  Data 


Figure  4.  Imaginary  (dielectric 

constant)  vs.  volumetric 
moisture  content, 

100  MHz,  20  degC 
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Figure  5.  Loss  tangent  vs. 

volumetric  moisture 
content,  100  MHz, 
20  deg  C 


All  ERDC  Samples 
100  MHz  Data 


volumetric  moisture  -  % 


All  ERDC  Samples 
100  MHz  Data 


Figure  6.  Apparent  conductivity 
vs.  volumetric  moisture 
content,  100  MHz, 

20  deg  C 


volumetric  moisture  -  % 
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Figure  7.  EM  power  attenuation 
vs.  volumetric  moisture 
content,  100  MHz, 

20  deg  C 


All  ERDC  Samples 
100  MHz  Data 


All  ERDC  Samples 
100  MHz  Data 


Figure  8. 


Normalized  phase 
velocity  vs.  volumetric 
moisture  content, 

100  MHz,  20  degC 
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dielectric  permittivity.  Since  permittivity  is  highly  correlated  with  moisture 
content,  then  it  should  be  no  surprise  that  phase  velocity  is  also  highly  correlated. 

In  summary,  single-frequency  soil  electrical  property  data  (100  MHz)  show 
that  dielectric  permittivity  and  phase  velocity  are  strong  functions  of  soil 
volumetric  moisture  content.  All  other  data,  including  EM  power  attenuation,  are 
not  strongly  correlated  with  moisture  content.  Further  analyses  in  Chapter  3  will 
focus  on  other  properties  of  the  soils,  including  measurement  frequency,  soil 
texture,  bulk  density,  and  physical  temperature,  in  an  effort  to  produce  models 
that  can  predict  soil  EM  attenuation  from  relatively  simple  field  measurable 
parameters. 


Soil  Classification  Schemes 

The  goal  of  this  study  is  to  derive  EM  attenuation  prediction  models  for  soils 
using  the  least  number  of  parameters  possible  that  can  be  easily  obtained  by  an 
experienced  soil  scientist  or  soils  engineer  working  in  the  field.  It  is  believed 
that  one  of  those  parameters  would  be  a  soil  type  classification  label.  Two 
different  soil  classification  systems  were  used  in  this  study  to  label  as  many  of 
the  ERDC-tested  soil  samples  as  possible.  One  of  them  came  from  the 
agricultural  community,  and  the  other  grew  out  of  the  need  to  label  the 
engineering  properties  of  soils. 


Soil  texture  triangle 

The  U.S.  Department  of  Agriculture  (USDA)  soil  classification  system  is  an 
adaptation  of  a  late- 19*  century  Russian  system  that  permitted  the  study  of  soils 
with  the  same  agricultural  characteristics  (Portland  Cement  Association  1973). 
In  other  words,  agricultural  experts  desired  a  way  to  classify  soils  so  that  soils  in 
various  parts  of  the  world  that  possessed  the  same  drainage  characteristics  and 
that  came  from  the  same  topography  would  grow  the  same  kinds  of  crops  when 
exposed  to  similar  weather  conditions.  This  classification  scheme  leads  to 
labeling  a  soil  with  a  great  soil  group  name  such  as  alluvial  soils  or  laterite  soils. 
The  great  soil  group  name  reflects  the  soil  maturity  and  the  climate,  vegetation, 
and  topography  of  its  location.  Each  group  is  further  divided  into  series.  All 
soils  within  a  series  are  similar  in  all  respects  except  the  texture.  Texture  can  be 
defined  as  the  relative  distribution  of  particles  sizes  within  a  soil  sample.  It  is 
this  final  textural  classification  that  was  used  to  label  the  ERDC  soil  samples. 

In  graphical  form.  Figure  9  defines  the  12  basic  textural  groups  used  by  the 
USDA  Natural  Resources  Conservation  Service.  To  properly  use  this  chart,  one 
first  needs  to  calculate  the  relative  amounts  of  specified  particle-size  bands. 
Using  only  that  portion  of  a  soil  that  will  pass  through  a  No.  10  sieve  (2  mm), 
a  particle-size  distribution  curve  can  be  generated  in  a  laboratory.  The  portion 
called  sand  consists  of  particles  between  0.05  and  2.0  mm  in  size;  silt  falls 
between  0.002  and  0.05  mm;  and  clay  consists  of  everything  smaller  than 
0.002  mm  in  size.  Each  point  on  the  textural  triangle  in  Figure  9  represents  a 
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Figure  9.  U.S.  Department  of 

Agriculture  soil  texture 
triangle 


unique  combination  of  sand  and  clay  content.  For  a  given  textural  class,  all 
combinations  of  sand  and  clay  content  are  bound  by  a  polygon  that  bears  the 
name  of  the  class.  For  example,  a  soil  possessing  30  percent  sand  and  35  percent 
clay-sized  particles  would  be  classified  as  a  “clay  loam.” 


Engineering  classification 

In  1942  the  U.S.  Army  Corps  of  Engineers  adopted  the  “Airfield 
Classification”  system  for  soils  that  had  been  developed  by  Arthur  Casagrande,  a 
professor  of  civil  engineering  at  Harvard  University.  With  minor  alterations  it 
became  what  today  is  known  as  the  Unified  Soil  Classification  System  (USCS). 
Its  purpose  was  to  classify  soils  in  a  way  that  would  indicate  how  the  soil  would 
behave  as  an  engineering  construction  material  (U.S.  Army  Engineer  Waterways 
Experiment  Station  1960). 

Eike  the  USDA  soil  texture  triangle,  the  USCS  is  based,  in  part,  on  the  shape 
of  the  soil  particle-size  distribution  curve.  However,  the  USCS  technique 
considers  all  sizes  of  material  within  the  soil.  Furthermore,  it  labels  sands  as 
having  grain  sizes  between  0.074  mm  (No.  200  sieve)  and  4.76  mm  (No.  4 
sieve),  while  silts  and  clays  are  those  particles  smaller  than  0.074  mm.  The 
USCS  goes  on  to  include  in  its  classification  procedure  two  other  laboratory- 
measured  parameters  called  the  liquid  limit  and  the  plastic  limit  (which  are 
related  by  a  term  called  the  plasticity  index),  which  are  collectively  known  as  the 
Atterberg  limits.  The  liquid  limit  is  the  weight-based  moisture  content  above 
which  the  soil  flows  as  a  viscous  liquid,  and  below  which  it  is  plastic.  The 
plastic  limit  is  the  moisture  content  at  which  the  soil  will  start  to  crumble  when 


14 


Chapter  2  ERDC  Soil  Data 


rolled  into  a  thread  under  the  palm  of  the  hand.  The  plasticity  index  is  just  the 
difference  between  the  liquid  limit  and  the  plastic  limit 

Classification  by  the  USCS  involves  following  a  decision  flowchart  that 
combines  key  particle-size  distribution  curves  values  with  the  soil  plasticity 
parameters  defined  above  to  arrive  at  one  of  15  different  classification  labels. 
Each  label  has  associated  with  it  a  word  description  of  the  soil.  For  example,  an 
OL  soil  is  described  as  “organic  silts  and  organic  silty  clays  of  low  plasticity.” 
The  USCS  silt  covers  a  different  particle-size  band  than  does  its  USDA 
counterpart. 

USDA  and  USCS  classification  labels  were  applied  to  all  of  the  ERDC  soil 
samples  for  which  such  information  was  available.  Of  the  1,080  data  samples 
shown  on  the  first  few  charts,  approximately  375  of  those  received  USDA 
classification  labels.  Approximately  460  could  be  classified  according  to  the 
USCS.  These  will  constitute  the  data  that  are  analyzed  in  the  following  chapter. 


ERDC  soils  mapped  onto  soil  triangle 

An  obvious  question  that  needs  to  be  answered  is  whether  or  not  USDA  and 
USCS  classifications  are  comparable  in  some  straightforward  way.  To  answer 
this  question,  all  of  the  soil  samples  tested  at  ERDC  that  could  be  classified  by 
both  methods  as  well  as  a  few  soils  whose  classifications  were  found  in  other 
reports  were  mapped  onto  the  USDA  soil  texture  triangle  shown  in  Figure  10.  If 
there  was  a  strong  correlation  between  the  two  classification  schemes,  then  each 
USCS  label  would  be  concentrated  in  only  one  or  two  of  the  USDA  triangle 
polygons.  Clearly,  this  is  not  the  case.  One  has  to  conclude  that  there  is  no 
simple  way  to  transfer  USDA  classifications  to  the  USCS  and  vice  versa.  Each 
scheme  will  need  to  be  examined  independently  for  its  ability  to  produce  an  EM 
attenuation  prediction  model. 
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•  CH.  high  plasticity  clay 

•  CL.  low  plasticity  clay 
▲  MH.  high  plasticity  silt 
A  ML.  low  piasticity  silt 

■  SC.  clayey  sand 

■  SM.  silty  sand 

•  SP.SW  sand 

•  GC ,  clayey  gravel 
GM .  silty  gravel 

•  GP.GW  gravel 
OH  organic  clay 
OL  organic  silt 


%  Silt 


%  sand 


Figure  10.  Mapping  of  USCS  ciassifications  onto  the  USDA  soii  triangie 
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3  Analysis  of  Attenuation 
Data 


Regression  Analyses 

As  noted  earlier,  the  ultimate  goal  of  this  study  was  to  develop  useful  EM 
power  attenuation  predietion  models  whose  input  would  eonsist  of  parameters 
easily  measured  in  the  field.  If  laboratory  measurements  were  required  to 
provide  model  input,  then  direet  attenuation  measurements  in  a  laboratory  eould 
be  made  just  as  easily,  and  predietive  models  would  not  be  neeessary. 

Analysis  of  the  ERDC  database  began  with  blind  linear  regressions  using 
SPSS  9. 0,  a  eomprehensive  eomputer  software  paekage  for  analyzing  data.  EM 
power  attenuation  was  taken  to  be  the  dependent  variable,  while  the  independent 
variables  were  ehosen  to  be  frequeney,  volumetrie  moisture  content,  and  sample 
bulk  density.  One  set  of  regressions  were  performed  for  the  USDA  soil 
classifications,  and  another  set  were  performed  for  the  USCS  soil  classifications. 
In  some  cases,  SPSS  could  not  produce  a  model;  in  others  the  fits  were  of  very 
poor  quality  as  evidenced  in  low  values  (R^  is  an  estimate  of  the  variance  of  a 
model  fit  from  actual  data  and  is  often  referred  to  as  the  residual  sum  of  squares 
or  the  error  sum  of  squares). 

It  was  quickly  established  that  fits  needed  to  be  attempted  at  fixed  values  of 
frequency.  Models  were  generated  at  frequencies  of  10  MHz,  100  MHz,  500 
MHz,  1,000  MHz,  5,000  MHz,  and  10,000  MHz  for  each  set  of  soil 
classifications.  Once  again,  results  were  very  poor.  Showing  tables  of  these 
results  would  be  meaningless  at  this  point  in  light  of  the  arguments  that  are  going 
to  be  made  in  the  following  sections. 

Clear  evidence  of  the  failure  to  produce  accurate  predictive  models  can  be 
seen  in  Figures  11-14.  These  are  charts  of  EM  power  attenuation  plotted  against 
volumetric  moisture  content  at  fixed  frequencies.  The  sample  temperature  was 
always  20  deg  C,  and  the  data  points  represent  the  entire  range  of  bulk  densities 
measured  at  ERDC.  Figures  1 1  and  12  are  for  USDA  classifycations  and  for  100 
MHz  and  1,000  MHz,  respectively.  Figures  13  and  14  represent  the  same 
parameters  but  for  USCS  classifications. 
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Figure  1 1 .  EM  power  attenuation  vs.  volumetric  moisture  at  1 00  MHz  for  USDA  classifications 
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Figure  1 2.  EM  power  attenuation  vs.  volumetric  moisture  at  1 ,000  MHz  for  USDA  classifications 
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Figure  1 3.  EM  power  attenuation  vs.  voiumetric  moisture  at  1 00  MHz  for  USCS  ciassifications 
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Figure  14.  EM  power  attenuation  vs.  voiumetric  moisture  at  1 ,000  MHz  for  USCS  ciassifications 
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If  there  was  a  possibility  of  produeing  EM  power  attenuation  predietion 
models  at  fixed  frequeneies  and  for  known  soil  types,  then  the  data  shown  on 
these  figures  would  be  elustered,  or  banded,  by  soil  elassifieation.  It  is  obvious 
from  examining  these  figures  that  this  is  not  the  ease.  All  of  the  soil 
elassifieation  symbols  are  mixed  throughout  eaeh  ehart. 

Another  goal  of  these  analyses  was  to  reduee  the  number  of  soil 
elassifieations  so  that  the  engineer  or  soil  seientist  in  the  field  would  not  have  to 
do  a  preeise  elassifieation.  For  example,  the  hope  was  that  sand,  loamy  sand,  and 
sandy  loam  USDA  elassifieations  eould  be  eonsolidated  into  one  elassifieation 
ealled  sandy  soil.  Loam,  silt  loam,  and  silt  might  eondense  into  a  silty  soil 
elassifieation.  All  of  the  other  soil  types  would  then  be  elayey  soils.  Needless  to 
say,  if  12  elassifieations  failed  to  produee  a  useful  predietion  model,  it  should  be 
no  surprise  that  three  elassifieations  would  also  fail.  The  same  kind  of  exereise 
was  performed  on  the  USCS  elassifieations  with  the  same  results. 

Although  soil  eonduetivity  predietive  models  have  been  reported  in  the 
literature  and  diseussed  in  some  detail  in  the  first  ehapter  of  this  report,  one  must 
keep  in  mind  that  those  models  were  generated  for  individual  soil  samples,  not 
for  a  eolleetion  of  samples  that  shared  a  eommon  elassifieation.  The  SPSS 
studies  of  the  ERDC  database  of  soil  eleetrieal  properties  elearly  established  that 
aeeurate  EM  power  attenuation  predietion  models  for  broad  elasses  of  soils  is  not 
aehievable. 

At  this  point  the  deeision  was  made  to  reexamine  the  ERDC  data  from  a 
qualitative  perspeetive  with  the  hope  of  establishing  some  sort  of  guidanee  that 
eould  still  be  useful  to  field  personnel.  These  exereises  follow. 


Bulk  Density  Effects 

Is  it  possible  that  some  elustering  eould  be  aehieved  by  dividing  the  data  into 
several  bands  of  bulk  density  values?  Figures  15  and  16  represent  plots  of  EM 
power  attenuation  vs.  volumetrie  moisture  eontent  for  all  of  the  ERDC  data  at 
two  different  frequeneies  but  where  the  data  have  been  broken  into  three  bulk 
density  bands;  namely,  values  less  than  1.4  g/ee,  values  between  1.4  and  1.7  g/ee, 
and  values  greater  than  1.7  g/ee.  Onee  again,  if  bulk  density  was  to  be  a 
eontrolling  faetor  in  predieting  EM  power  attenuation,  then  the  data  on  these 
eharts  should  be  elustered,  and  they  are  not.  All  that  one  ean  say  is  that  higher 
frequeneies  result  in  higher  attenuation  values. 


Temperature  Effects 

A  small  amount  of  ERDC  data  exists  for  temperatures  other  than  20  deg  C. 
While  some  of  those  results  have  been  published  previously  (Curtis  1993),  it  is 
probably  worth  presenting  all  of  the  ERDC  variable  temperature  data  here  to 
show  whether  or  not  soil  temperature  might  provide  some  guidanee  to  the  field 
engineer/seientist.  Unfortunately  the  results  are,  onee  again,  ineonelusive. 
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Figure  1 5.  Bulk  density  effects  on 
EM  power  attenuation  at 
100  MHz  and  20  deg  C 
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Figure  16.  Bulk  density  effects  on 
EM  power  attenuation 
at  1,000  MHz  and 
20  deg  C 
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Figures  17  and  18  contain  plots  of  EM  power  attenuation  vs.  volumetric 
moisture  content  for  those  few  soil  samples  tested  at  several  temperatures  and  at 
frequencies  of  1 00  MFIz  and  1 ,000  MFIz,  respectively.  Close  examination  of 
Figure  17  (100  MFIz)  leads  one  to  the  conclusion  that  increasing  soil  temperature 
results  in  increasing  power  attenuation.  However,  1 ,000  MHz  data  shown  in 
Figure  18  produces  the  opposite  conclusion.  The  most  conservative  guidance, 
then,  is  simply  that  increasing  soil  temperature  may  result  in  less  power 
attenuation.  Assigning  numbers  to  the  temperature  effect  is  futile  in  light  of  the 
failure  to  predict  attenuation  for  a  given  soil  type  and  system  frequency. 


Frequency  Effects 

Given  that  regression  analyses  failed  to  produce  useful  EM  power 
attenuation  predictive  models,  and  qualitatively  determining  that  neither  soil 
density  nor  soil  temperature  had  a  strong  and  predictable  impact,  all  that  remains 
is  to  examine  the  effects  of  measurement  frequency  on  power  attenuation. 

Perhaps  there  will  be  trends  in  the  data  that  may  be  useful  for  field  applications. 

Figures  19-24  contain  charts  of  EM  power  attenuation  plotted  against 
volumetric  moisture  content  for  several  fixed  measurement  frequencies:  10  MHz, 
100  MHz,  500  MHz,  1,000  MHz,  5,000  MHz,  and  10,000  MHz.  These  charts 
represent  all  of  the  ERDC  data  collected  at  those  frequencies  and  at  a  temperature 
of  20  deg  C.  At  each  frequency  one  chart  is  plotted  on  a  linear  scale,  while  the 
following  chart  is  plotted  on  a  semi-logarithmic  scale.  The  reason  for  plotting 
data  on  a  logarithmic  scale  is  to  provide  a  better  visualization  of  the  spread  of  the 
data.  Note  also  that  the  vertical  scales  have  been  changed  on  each  chart  to  best 
view  the  data. 

Consider  first,  all  of  the  ERDC  data  collected  at  1 0  MHz  and  shown  on 
Figure  19a.  About  590  data  points  are  plotted  on  Figure  19a.  Included  on  that 
chart  is  a  line  representing  a  linear  least  squares  fit  to  the  data.  Note  how  poor 
the  fit  is  (R^  =  0.1 12).  Clearly  the  linear  fit  could  not  be  used  as  a  predictive 
model,  because  it  would  seriously  underpredict  values  of  attenuation.  For 
example,  at  20  percent  volumetric  moisture,  the  model  would  predict  about  5 
dB/m  of  attenuation.  However,  it  is  clear  from  the  data  that  a  large  number  of 
soil  samples  demonstrated  far  less  attenuation. 

On  a  logarithmic  scale  (Figure  19b),  one  can  readily  see  a  lower  bound  on 
the  data.  This  is  potentially  a  useful  tool  in  the  sense  that  it  gives  field  personnel 
a  conservative  lower  estimate  on  power  attenuation.  For  the  case  of  a  sensor 
operating  at  1 0  MHz,  that  conservative  estimate  for  EM  power  attenuation  might 
be  roughly  1  dB/m,  which  is  not  a  significant  amount  of  attenuation.  Experience 
has  shown  that  low-frequency  systems  can  penetrate  much  deeper  into  soils  than 
do  high-frequency  systems.  These  attenuation  numbers  support  that  observation. 
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Figure  17.  Temperature  effects  on 
EM  power  attenuation  at 
100  MHz  and  20  deg  C 
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Figure  18.  Temperature  effects  on 
EM  power  attenuation 
at  1,000  MHz  and 
20  deg  C 
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Power  Attenuation  vs  Moisture  Content 
All  ERDC  Samples 
10  MHz  Data 


Figure  19a.  EM  power  attenuation 

vs.  volumetric  moisture  at 
10  MHz  and  20  deg  C 
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Figure  19b. 


EM  power  attenuation 
vs.  volumetric 
moisture  at  10  MHz 
and  20  deg  C  (log 
scale) 
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Power  Attenuation  vs  Moisture  Content 
All  ERDC  Samples 
100  MHz  Data 


Figure  20a.  EM  power  attenuation 

vs.  volumetric  moisture  at 
100  MHz  and  20  deg  C 
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Figure  20b. 


EM  power  attenuation 
vs.  volumetric 
moisture  at  100  MHz 
and  20  deg  C  (log 
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Power  Attenuation  vs  Moisture  Content 
All  ERDC  Samples 
500  MHz  Data 


Figure  21  a.  EM  power  attenuation  vs. 

volumetric  moisture  at  500 
MHz  and 
20  deg  C 
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Figure  21b. 


EM  power  attenuation 
vs.  volumetric 
moisture  at  500  MHz 
and  20  deg  C 
(log  scale) 
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Power  Attenuation  vs  Moisture  Content 
All  ERDC  Samples 
1,000  MHz  Data 


Figure  22a.  EM  power  attenuation  vs. 

voiumetric  moisture  at 
1,000  MHz  and 
20  deg  C 
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Figure  22b. 


EM  power  attenuation 
vs.  voiumetric 
moisture  at  1,000  MHz 
and  20  deg  C 
(iog  scaie) 
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Power  Attenuation  vs  Moisture  Content 
All  ERDC  Samples 
5,000  MHz  Data 


Figure  23a.  EM  power  attenuation  vs. 

volumetric  moisture  at 
5,000  MHz  and 
20  deg  C 


Figure  23b. 


EM  power  attenuation 
vs.  volumetric 
moisture  at  5,000  MHz 
and  20  deg  C 
(log  scale) 
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Figure  24a.  EM  power  attenuation 
vs.  volumetric  moisture 
at  10,000  MHz  and 
20  deg  C 


Figure  24b. 


EM  power  attenuation 
vs.  volumetric 
moisture  at  10,000 
MHz  and  20  deg  C 
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Figures  20a  and  20b  show  attenuation  data  at  100  MFlz  (about  1,080  data 
points).  As  before,  the  linear  regression  model  is  not  good  (R^  =  0.2033).  Note 
also  that  the  average  power  attenuation  at  1 00  MFIz  is  more  than  twiee  the 
attenuation  at  10  MFIz.  A  elose  examination  of  the  logarithmie  ehart  (Figure  20b) 
reveals  at  least  two  interesting  eharaeteristies  of  the  attenuation  data.  One  is  that 
at  low  moisture  eonditions  (less  than  5  pereent  volumetrie  moisture),  the  soils 
tested  exhibited  a  huge  variation  in  power  attenuation;  however,  if  one  needs  to 
be  eonservative,  the  lower  limit  of  that  attenuation  is  eertainly  quite  low.  The 
seeond  observation  is  that  if  one  wanted  to  plaee  a  fairly  eonservative  estimate  on 
power  attenuation  at  normal  moisture  eontents  and  for  mueh  wetter  soils,  that 
number  eould  easily  be  ehosen  at  a  value  of  say,  2  or  3  dB/m. 

Figures  21a  and  21b  show  500  MFIz  of  data  (again  representing  about 
1,080  data  points).  The  linear  regression  model  is  no  better  at  500  MFIz  than  it 
was  at  100  MFIz;  however,  the  average  attenuation  has  nearly  doubled  for  the 
higher  frequeney  data.  Other  than  for  the  lowest  values  of  moisture  eontent,  one 
eould  argue  that  a  reasonably  eonservative  estimate  of  power  attenuation  might 
be  on  the  order  of  5  dB/m  at  500  MFIz. 

The  linear  regression  model  for  1,000  MFIz  data  (Figures  22a  and  22b)  is  a 
slight  improvement  over  the  500  MFIz  model,  but  it  is  still  not  a  useful  predietive 
tool.  A  eonservative  attenuation  value  of  1 0  dB/m  or  more  is  possible  for  normal 
to  high  moisture  eontents.  Figures  22a  and  22b  also  eontain  about  1,080  data 
points. 

It  is  not  until  one  examines  the  higher  frequeney  data  at  5,000  MFIz  (Figures 
23  a  and  23b)  that  one  sees  a  reasonably  good  linear  regression  model  (R  = 
0.9299).  There  are,  however,  a  number  of  eaveats  that  must  be  stipulated  along 
with  that  observation.  One  is  that  far  fewer  data  points  existed  for  the  5,000 
MHz  frequeney  (approximately  150).  In  faet,  these  data  represent  soils  for  one 
major  field  sampling  eolleetion,  along  with  about  a  dozen  “laboratory”  soils 
whose  partiele-size  distribution  eurves  were  highly  variable,  but  whose  ehemistry 
differenees  might  not  have  been  extreme. 

The  same  arguments  hold  for  the  10,000  MHz  data  (Figures  24a  and  24b), 
whieh  represent  only  about  65  data  points  and  are  limited  mostly  to  the 
“laboratory”  soils.  In  other  words,  all  that  ean  be  said  about  the  highest 
frequeney  data  is  that  average  attenuation  values  eontinue  to  rise  as  the  frequeney 
goes  up.  Care  should  be  taken  in  extrapolating  those  results  to  dry  soils.  All  of 
the  data  elearly  show  that  for  the  driest  of  soils,  one  ean  expeet  signifieant 
penetration  of  the  EM  signal. 
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4  Summary  and 
Recommendations 


The  goal  of  this  study  was  to  determine  if  EM  power  attenuation  predietion 
models  eould  be  developed  for  soils  that  required  only  a  minimum  of  easily 
obtained  input  parameters.  Field  erews  using  high-power  eleetromagnetie 
deviees  to  deteet  and/or  diseriminate  buried  unexploded  ordnanee  (UXO)  need  to 
know  whether  or  not  the  power  levels  of  their  equipment  are  suffieient  to 
detonate  UXO  fuzes.  Knowing  the  power  level  of  their  transmitters  and  the 
attenuation  faetor  for  the  soil  they  are  penetrating,  they  eould  prediet  power 
levels  at  the  antieipated  depths  of  UXO. 

Faetors  that  affeet  EM  attenuation  inelude  sensor  operating  frequeney,  soil 
ehemistry,  moisture  eontent,  soil  temperature,  and  bulk  density.  Based  on  the 
author’s  experienee  in  measuring  soil  eleetrieal  properties,  the  original  hope  for 
this  study  was  that  useful  models  eould  be  exeeuted  using  soil  type,  frequeney, 
and  moisture  eontent  as  a  minimum  set  of  input  parameters. 

Sinee  the  early  1990s  eomplex  dieleetrie  properties  of  more  than  300  distinet 
soil  samples  have  been  measured  at  ERDC  in  Vieksburg,  MS.  A  dozen,  or  so,  of 
these  soils  were  obtained  from  the  ERDC  Geoteehnieal  Eaboratory;  and,  while 
they  eannot  be  elassified  as  artifieial  soils,  they  were  eertainly  proeessed  to 
remove  organie  material  and  gravel-sized  partieles.  Another  group  of  five  soils 
eame  from  the  U.S.  Natural  Resourees  Conservation  Serviee  and  ean  be  referred 
to  as  “benehmark”  soils.  The  remainder  of  the  soils  were  all  obtained  from 
numerous  field  sites  in  this  eountry  as  well  as  from  some  foreign  sites.  A  variety 
of  frequeney  bands  were  required  for  different  eustomers.  Some  data  were 
eolleeted  at  several  sample  temperatures;  most  were  eolleeted  at  20  deg  C.  The 
dieleetrie  properties  of  eaeh  sample  were  measured  for  at  least  three  different 
volumetrie  moisture  eontents.  Sample  density  eould  not  be  eontrolled  in  these 
experiments;  therefore,  a  broad  range  of  random  sample  densities  were  aehieved. 
In  summary,  at  seleeted  frequeneies,  more  than  1 ,000  EM  power  attenuation  data 
points  were  available  for  modeling  studies. 

Not  all  of  the  soils  tested  at  ERDC  were  elassified  at  the  time  that  their 
dieleetrie  properties  were  being  measured.  However,  the  author  was  able  to 
reeover  one  of  two  different  elassifieations  for  many  of  the  samples,  and  both 
elassifieations  for  a  number  of  samples.  One  elassifieation  seheme  is  that  used 
by  the  U.S.  Department  of  Agrieulture.  Naturally  its  objeetive  is  to  be  able  to 
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compare  soils  from  different  parts  of  the  world  in  terms  of  their  ability  to  sustain 
crop  growth.  The  USDA  classification  scheme  can  be  visualized  as  a  soil  texture 
triangle,  within  which  a  soil  is  classified  according  to  its  relative  percentages  of 
sand  and  clay-sized  particles.  The  second  classification  scheme  is  one  used  by 
the  U.S.  Army  Corps  of  Engineers  and  referred  to  as  the  Unified  Soil 
Classification  System.  USCS  soils  are  labeled  according  to  their  performance  as 
load-bearing  construction  materials. 

The  initial  effort  at  developing  EM  power  attenuation  prediction  models 
involved  numerous  attempts  to  perform  multiple  parameter  linear  regression 
analyses  of  the  ERDC  data.  These  efforts,  which  met  with  no  success,  focused 
on  trying  to  develop  predictive  models  for  each  soil  classification  at  set 
frequencies,  a  single  sample  temperature,  and  utilizing  moisture  content  and  bulk 
density  as  input  parameters.  Data  charts  were  shown  in  Chapter  3  that  clearly 
demonstrated  that  soil  type  could  not  be  used  to  consistently  or  definitively 
predict  EM  power  attenuation.  Similar  charts  revealed  that  neither  sample 
density  nor  sample  temperature  provided  any  data  clustering  when  all  of  the  soil 
types  were  displayed  on  a  single  chart. 

Following  the  failure  of  the  statistical  analyses  to  generate  the  desired 
predictive  models,  the  total  ERDC  data  set  was  revisited  at  each  of  five 
frequencies,  and  a  simple  linear  regression  was  performed.  High-frequency  data 
(>1  GHz)  showed  reasonably  good  model  fits,  but  probably  only  because  the 
soils  tested  at  those  frequencies  were  of  the  processed  variety.  It  was  concluded, 
however,  that  the  data  charts  for  all  of  the  frequencies  could  be  used  to  provide  a 
conservative  lower  bound  on  power  attenuation.  In  light  of  that  observation,  the 
following  table  has  been  generated  from  the  data  charts  shown  in  Figures  19-24. 
For  each  of  the  six  selected  frequencies,  the  linear  regression  model  was  used  to 
identify  an  average  attenuation  at  1 0  percent  volumetric  soil  moisture  content. 
The  figure  of  1 0  percent  was  selected  by  the  author  as  a  typical  in  situ  moisture 
value  for  near-surface  soils  at  most  temperate  locales.  Upon  examining  the 
spread  of  the  data  at  each  frequency  and  at  the  1 0  percent  moisture  value,  a 
column  of  conservative  attenuation  values  was  generated  and  included  in  the 
table.  These  are  values  that  would  provide  some  degree  of  safety  for  field  crews 
operating  high-power  sensing  equipment  in  the  presence  of  UXO. 


Table  1 

Frequency-Dependent  EM  Power  Attenuation  Values 

Average  Attenuation  at  10% 

Conservative  Attenuation 
Estimate  at  10%  Volumetric 

Frequency  (MHz) 

Volumetric  Moisture  (dB/m) 

Moisture  (dB/m) 

10 

3 

1 

100 

6 

2 

500 

12 

4 

1,000 

17 

10 

5,000 

122 

100 

10,000 

422 

300 
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While  the  information  in  Table  1  may  be  useful  guidanee  in  a  general  sense, 
it  does  not  substitute  for  laboratory  evaluation  of  test  site  soil  properties.  The 
only  definitive  approaeh  to  providing  a  safe  work  environment  is  to  eolleet 
numerous  near-surfaee  and  below-surfaee  samples  and  transport  them  to  a 
laboratory  where  EM  power  attenuation  measurements  ean  be  made  at  the 
appropriate  frequeney. 
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Appendix  A 

Soil  Dielectric  Property  Data 
Collection  and  Analysis 
Procedures 


Complex  Dielectric  Property  Measurements 

From  about  1990  until  the  present,  the  eomplex  dieleetrie  eonstant  of  soils 
has  been  measured  by  the  author  at  the  U.S.  Army  Engineer  Researeh  and 
Development  Center,  Vieksburg  Site,  utilizing  a  Flewlett-Paekard  8753E 
Network  Analyzer  System  (and  an  earlier  FIP8510C)  with  a  eoaxial  sample 
holder  (Curtis  2001b).'  The  idealized  geometry  used  for  these  measurements  is  a 
sample  of  length,  d,  in  a  eoaxial  holder.  The  eomplex  relative  dieleetrie  eonstant 
and  eomplex  relative  magnetie  permeability  of  a  sample  are  labeled  s  and  /^, 
respeetively.  Complex  S-parameter  data  are  eolleeted  at  seleeted  ffequeneies 
and  stored  on  a  laboratory  data  eontroller.  The  parameter,  Su,  is  the  eomplex 
ratio  of  refleeted  voltage  to  ineident  voltage  at  one  end  of  the  sample  (with  the 
other  end  eleetronieally  shorted),  while  the  parameter,  S21,  is  the  eomplex  ratio  of 
voltage  measured  on  the  other  end  of  the  sample  to  the  same  ineident  voltage  (an 
idealized  transmission  eoeffieient). 

The  use  of  S-parameters  as  measurement  tools  for  eleetronie  deviees  was 
introdueed  in  the  1960s  (Flewlett  Paekard).  In  1970,  a  paper  was  published  that 
deseribed  the  use  of  eoaxial  samples  and  S-parameter  measurements  to  determine 
the  eomplex  dieleetrie  and  magnetie  properties  of  a  material  (Nieolson  and  Ross 
1970).  The  model  ean  be  summarized: 


and 


’  References  cited  in  this  appendix  are  located  at  the  end  of  the  main  text. 
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c  is  the  speed  of  light  in  a  vaeuum,  and  m  is  the  angular  ffequeney  of  the  signal. 
Implied  in  this  model  is  the  assumption  that  the  sample  holder  has  no  intrinsie 
magnetie  properties.  This  is  undoubtedly  not  true  for  the  sample  holders  used  in 
this  projeet.  In  faet,  the  eenter  eonduetor  for  the  sample  holders  are  made  of  a 
steel  alloy,  and  the  eonneetors  at  the  ends  of  the  sample  holders  are  not  preeision, 
low-loss  deviees.  Nevertheless,  measurements  of  empty  sample  holders  (air)  and 
holders  filled  with  distilled,  deionized  water  give  very  aeeeptable  results.  If  one 
is  able  to  eonfidently  assume  that  the  material  being  measured  is  nonmagnetie, 
then  the  eomplex  dieleetrie  eonstant  ean  be  ealeulated  direetly  from  Equation  A4. 


Calculation  of  Additional  Parameters 

While  the  eomplex  dieleetrie  eonstant  and  eomplex  magnetie  permeability 
are,  by  themselves,  meaningful  properties  of  any  material,  it  is  possible  to 
generate  from  them  other  useful  parameters  that  are  more  readily  utilized  by 
engineers  and  seientists.  For  example,  in  the  design  of  ground-penetrating  radar 
systems,  what  may  be  most  useful  to  the  engineer  is  how  rapidly  the  soil 
attenuates  the  signal.  Attenuation,  eoupled  with  the  sensitivity  of  the  reeeiver, 
determines  the  effeetive  depth  of  penetration  for  the  radar  in  a  partieular  soil. 


A2 
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Similarly,  a  clear  understanding  of  how  the  speed  of  an  electromagnetic  wave 
varies  with  depth  and  lateral  position  at  any  given  test  site  determines  how  well 
the  radar  data  can  be  processed  to  produce  an  accurate  assessment  of  subsurface 
conditions.  The  following  paragraphs  outline  the  generation  of  several  useful 
parameters  from  complex  dielectric  properties,  including  power  attenuation, 
phase  velocity,  and  an  equivalent  electrical  conductivity. 

Assuming  plane  harmonic  wave  propagation  in  a  lossy,  unbounded  medium, 
the  wave  amplitude  function  may  be  written: 

^i(kx  -  cot) 

where 


A:  =  P  +  ic(=  (i)N/c 

P 

a 

00 

N 

c 


=  complex  propagation  constant, 
=  phase  constant, 

=  amplitude  attenuation  factor, 

=  radial  frequency, 

=  complex  index  of  refraction, 

=  velocity  of  light  in  a  vacuum. 


/■  =  symbol  designating  an  imaginary  quantity  =  , 

X  =  space  coordinate,  and 
t  =  time 

For  the  sake  of  simplicity,  let  the  medium  be  nonmagnetic.  Then, 


=  E  =  e'  ^  is" 


(AlO) 


where,  as  stated  above,  s  is  the  relative  complex  dielectric  constant.  The 
dielectric  constant,  along  with  the  electrical  conductivity  from  Ohm's  Law, 
represents  the  electrical  properties  of  the  medium.  Conductivity,  a ,  accounts  for 
current  flow  in  the  sample  because  of  free  charged  particle  motion;  whereas,  the 
real  part  of  the  dielectric  constant,  s' ,  accounts  for  Maxwell’s  “displacement 
currents”  because  of  the  electric  polarization  of  the  medium  and  to  the  time  rate- 
of-change  of  the  electric  field.  When  both  conduction  and  displacement  currents 
are  considered,  one  finds  two  terms  in  Ampere's  law  for  current  flow  that 
represent  losses  (or  a  shift  in  phase),  one  containing  the  electrical  conductivity 
and  one  containing  the  imaginary  part  of  the  dielectric  constant.  While  these  two 
terms  account  for  different  loss  mechanisms,  most  researchers  use  only  one  term 
or  the  other  to  identify  losses.  Many  users  prefer  to  deal  with  the  concept  of 
electrical  conductivity.  Flowever,  the  algorithm  used  in  this  laboratory  to 
measure  dielectric  properties  assumes  that  all  losses  are  contained  in  the  f  "term. 
Upon  observing  how  the  two  loss  terms  appear  in  Maxwell’s  equations,  one  can 
choose  to  use  them  interchangeably  through  the  relationship  (in  MKS  units). 
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(All) 


cr  =  a"  So  CO 

where  the  units  of  eonduetivity  are  mhos/meter  (or  siemens/meter)  and  So  is  the 
permittivity  of  free  spaee  (8.85x10'^^  farads/meter). 

Squaring  the  expression  for  the  eomplex  propagation  eonstant,  substituting 
the  expression  for  the  square  of  the  eomplex  index  of  refraetion,  and  equating 
real  and  imaginary  eomponents,  one  obtains  two  algebraie  equations  that  relate 
the  amplitude  attenuation  faetor  and  phase  eonstant  to  the  eomplex  dieleetrie 
eonstant: 


nl  2 

p  -  a 


(A12) 


and 
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(A13) 


Solving  these  equations  for  the  amplitude  attenuation  faetor  and  for  the  phase 
eonstant  results  in  the  following  expressions:’ 
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(A14) 


(A15) 


’  Had  magnetic  permeability  been  included,  then  the  phase  constant  would  have  been: 

(e'p'-£"p")  ±  +(£'P")"  +(£>')"  +(£>")" 

and  the  attenuation  coefficient  would  have  been: 

(B^  (  e'p"  +  e"p') 


2c"  L 
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2^c 
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The  f  %'ratio  is  also  referred  to  as  the  loss  tangent.  Keep  in  mind  that  some 
researehers  prefer  to  work  with  the  eleetrieal  eonduetivity  in  plaee  of  the 
dieleetrie  loss  term,  so  that  the  radieand  might  look  different  in  other  referenees 
(e.g.,  s'Vs'  may  be  replaeed  by  a/(s'soco)). 

Plane  waves  of  eonstant  phase  will  propagate  with  a  veloeity 

1 
2 

(A16) 


This  phase  veloeity  is  not  neeessarily  the  speed  with  whieh  the  energy  of  the 
wave  propagates  through  the  medium.  The  latter  is  referred  to  as  the  group 
veloeity  and  ean  be  ealeulated  as  the  rate  of  ehange  of  radial  frequeney  with 
respeet  to  the  phase  eonstant.  However,  as  long  as  the  phase  veloeity  is  relatively 
eonstant  over  the  range  of  ffequeneies  of  interest,  then  there  is  little  differenee 
between  phase  veloeity  and  group  veloeity. 

The  power  intensity  of  the  plane  eleetromagnetie  wave  deereases 
exponentially  with  depth  of  penetration  by  the  faetor,  e  ,  or,  in  one  unit  of 
distanee  traveled,  a  deerease  of  e  Power  attenuation  expressed  in  deeibels 
per  meter  ean  then  be  written  as: 
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